Astronomy & Astrophysics manuscript no. toloba'v2 


©ESO 2010 


December 3, 2010 





Formation and evolution of dwarf early-type galaxies in the Virgo 

cluster 

I. Internal kinematics 

E. Toloba 1 , A. Boselli 2 , A. J. Cenarro 3 , R. F. Peletier 4 , J. Gorgas 1 , A. Gil de Paz 1 , and J. C. Munoz-Mateos 15 

1 Departamento de Astrofisica y CC. de la Atmosfera, Universidad Complutense de Madrid, 28040, Madrid, Spain 

e-mail: etj@astrax.fis.ucm.es e-mail: jgorgas@fis.ucm.es e-mail: gildepaz@gmail.com e-mail: 
j cmunoz@astr ax . f i s . ucm . es 

2 Laboratoire d'Astrophysique de Marseille, UMR 61 10 CNRS, 38 rue F. Joliot-Curie, F-13388 Marseille, France 
e-mail: alessandro . boselli@oamp . fr 

3 Centro de Estudios de Ffsica del Cosmos de Aragon, E-44001, Teruel, Spain 
e-mail: cenarro@cefca.es 

4 Kapteyn Astronomical Institute, Rijksuniversiteit Groningen, Postbus 800, 9700 AV Groningen, the Netherlands 
e-mail: peletier@astro . rug . nl 

5 National Radio Astronomy Observatory, 520 Edgemont Road, Charlottesville, VA 22903-2475 

Received ; accepted 

ABSTRACT 

We present new medium resolution kinematic data for a sample of 21 dwarf early-type galaxies (dEs) mainly in the Virgo cluster, 
obtained with the WHT and INT telescopes at the Roque de los Muchachos Observatory (La Palma, Spain). These data are used 
to study the origin of the dwarf elliptical galaxy population inhabiting clusters. We confirm that dEs are not dark matter dominated 
galaxies, at least up to the half-light radius. We also find that the observed galaxies in the outer parts of the cluster are mostly 
rotationally supported systems with disky morphological shapes. Rotationally supported dEs have rotation curves similar to those of 
star forming galaxies of similar luminosity and follow the Tully-Fisher relation. This is expected if dE galaxies are the descendant of 
low luminosity star forming systems which recently entered the cluster environment and lost their gas due to a ram pressure stripping 
event, quenching their star formation activity and transforming into quiescent systems, but conserving their angular momentum. 

Key words. Galaxies: clusters: individual: Virgo Galaxies: dwarf Galaxies: elliptical and lenticular, cD Galaxies: kinematics and 
dynamics Galaxies: evolution Galaxies: dark matter 
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I~J I 1. Introduction pact and low su rface brightness dwarfs) (Bender et al.l (Il992l) : 

^ ■ iKormendv etal] (2009). Later on, it was shown that dEs were 

■ The processes involved in galaxy formation and evolution no longer small Es with simple, old and metal-poor stellar pop- 

; through cosmic time are still poorly understood. It is indeed still ulations, but much more complex objects exhibiting a wide 

O . unclear how matter assembled to form the present day galaxy range of stellar contents. For example, in the Virgo cluster, they 

^ ; population, whether it followed a passive evolution after the have stellar populations ranging from very young (around 1 Gyr 

>. ! collapse of the primordial density fluctuations (secular evolu- old) luminosity -weighted ages to as old as the oldest Es galax- 

tion), through a subsequent merging of growing structures (hi- ies (14 Gyr) ( Michiels en et alj|2008l) . Their proximity allowed 



X 



erarchical formation) or a combination of the two. A way of detailed studies of their structural properties which indicated 

^ , quantifying the relative role of these different mechanisms is that, behind their elliptical appearance, dEs show a great variety 

■ - ■ to stu dy dwarf galaxies, the most numerous objects in the uni- of underlying structures, like discs, spiral arms, irregular fea- 

verse dFerguson & Binggelil [ l994l) . Their importance resides in tures, etc, making t he m a v ery heterogeneous class of galaxies 

the fact that these low-luminosity systems are expected to be the (Lisker et al. 2006a b, 20Q3). 

building blocks of massive galaxies in lambda cold dark mat- These evidences indicate a complex formation process shap- 

terfACDM) hierarchical merging scenarios (e.g.,| White & Reesj ing me evo i ution of dEs in c i usters . Two main different pro- 

| 1978t | White & Frenk| | 1991D . cesses have been proposed in the literature: the first mecha- 

Among dwarf galaxies, quiescent dwarfs (which we here nism is based on the idea that dEs are formed through in- 

define to be all quiescent galaxies with Mg > -18, including ternal processes, like supernova feedback, where the interstel- 

both dwarf ellipticals and spheroidals, hereafter indicated as lar medium (ISM) of the progenitor star forming galaxy is 

dEs) are of particular intere st since they are the most numer- swept away by the kine tic pressure generated by supernovae 

ous population in clusters (Ferguson & Binggelil ll994h . These (lYoshii & Arimotoll 19871) . althou gh it seems highly unlikely i n 

objects were originally thought to be the low luminosity ex- dark-matter dominated systems (ISilich & Tenorio-T^ gle 2001); 

tension of giant ellipticals (Es). Since the 1990s it is known the second mechanism rests upon external processes induced by 

that dEs are composed of several families of objects (e.g. com- the interaction with the hostile environment in which dEs reside 
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dSandage et al.lll985b iBlanton et alll2005l) . In a dense environ- 
ment several mechanisms are affecting galaxies. This might hap- 
pen through interactions with the intergalactic medium (IGM), 
as ram-pressure stripping (Boselli et al. 2008a b), galaxy-galaxy 
interaction s (e.g., [B yi'd & Valtonen 1990) and galaxy harass- 
ment (e.g., iMoore et al.lll998t iMastropietro et alj|2005l) . It has 
been shown that all these interactions are able to reproduce 
some of the observational prope rties of local dEs in clu sters, 
like their structural para meters dLisker et al.l 2006bl 120071) or 
their stellar populations jGeha et al.ll2002 , 120031: van Zee et al.l 
I2004bt iMichielsen etal]|2008t iPaudel et al.ll2010l) , but none of 
them have been tested versus the kinematic properties. In this 
regard, whereas in the harassment scenario the system is rapidly 
heated, leading to an increase of the velocity dispersion of th e 
galaxy and a decrease of its rotation ( Mastropie tro et alj |2005). 
in a ra m-pressure stripping event the angular momentum is con- 
served dBoselli et alf 2008a ,b). 

With the aim of using kinematic data to constrain dwarf 
galaxy evolution, we recently started an ambitious observational 
program at the Observatory El Roque de los Muchachos (La 
Palma, Spain) to gather medium resolution spectroscopic data of 
dEs in the Virgo cluster. In this paper we present a detailed analy- 
sis of the internal kinematics focusing our attentio n into the most 
rapidly rotating systems. In Toloba et al. (2009) we have stud- 
ied the kinematics as a function of local environment inside the 
Virgo cluster. Several typical scaling relations of pressure sup- 
ported systems, such as the Fundamental Plane as well as their 
stellar population properties will be analysed in a forthcoming 
communication. 

This paper is structured as follows: in Sections 2, 3 and 4 
we describe the sample selection, the observations and the data 
reduction process. In Section 5 we report the kinematic measure- 
ments paying special attention to the method employed and the 
stellar templates used. We also describe the procedure followed 
to obtain the radial kinematic profiles (Section 5.1), the central 
velocity dispersion and the maximum rotational speed of the se- 
lected galaxies (Section 5.2), making comparisons with previous 
works (Section 5.3). Combined with photometric data (Section 
6), the present kinematic observations are used to study the prop- 
erties of rotationally supported systems in the framework of var- 
ious models of galaxy evolution (Sections 7, 8 and 9). 

2. The sample 

The sample analysed in this work is composed of galaxies with 
M r > -1 6 classified as dE or dSO in the Virgo Cluster Catalog 
(VCC)by Bing geli et al.l (1985). All galaxies have been selected 
to have SP SS imaging and to be within the GALEX MIS fields 
(Bos elli et al.l l2005). thus to have a measured UV magnitude or 
an upper limit. To these we added a few field quiescent dwarfs 
(originally used as fillers in our observing runs) useful for com- 
parison in a statistical study. Out of the 43 Virgo galaxies satis- 
fying these requirements in the VCC, 18 have been observed for 
this work. To make the observations accessible to 2. 5-4. 2m tele- 
scopes, we chose those objects with the highest surface bright- 
ness. 

The field sample consists of early-type dwarfs in low den- 
sity regions with magnitudes between -18.5 < M r < -14.5 and 
distances similar to Virgo (375 km s~'< v < 1875 km s , 5-25 
Mpc). Quiescent objects have been selected assuming the colour 
criterion FUV-NUV > 0.9 or u - g > 1.2 when UV detections 
were not available. We observed only 3 out of 10 field dEs candi- 
dates. To these 18 Virgo and 3 field dEs we added M32, selected 
to test the setup of the instruments. Thanks to the large amount 



of available data, M32 is also an ideal target for comparison with 
other data available in the literature. 

Although we can consider it representative of the bright end 
of the Virgo Cluster dE population, the observed sample is not 
complete in any sense. 



3. Observations 

The observing time that we obtained for this work was part of 
the International Time Program (ITP 2005-2007) at El Roque 
de los Muchachos Observatory. Here we focus on the medium 
resolution (R ^ 3800), long-slit spectroscopy carried out during 
three observing runs. In runs 1 and 3 (December 2005, February 
2007) we used the ISIS double-arm spectrograph at the 4.2m 
WHT, and in run 2 (January 2007) we used the IDS spectrograph 
at the INT (2.5m telescope). 

The advantage of ISIS over IDS is that it allows us to use 
a dichroic (5300 dichroic in our case) to split the light into two 
beams to observe simultaneously two wavelength ranges, one in 
the blue optical part of the spectra and another in the red. This 
technique allowed us to cover, in 3 settings in the first run, the 
full wavelength range from 3500 A to 8950 A, using a mirror 
to cover 5000-5600 A, the only range that we could not cover 
with this dichroic. In the third run we used 2 settings to cover 
the same wavelength range except for the dichroic gap. 

The wavelength range covered by the IDS was smaller 
(4600-5960 A), since detector and grating are the same as on the 
blue arm of ISIS, the data obtained had similar resolution. The 
spectral resolution (R ^ 3800) is high enough to obtain reliable 
kinematics for dwarf galaxies. 

All the details of the configurations used in each run are 
specified in Table Q] 

In Table [2] we list the observed sample. Colum n 5 presents 
the m orpholog ical type classifica tion according to Liske r et al.l 
(2006a!) and lLiskeretall d2006bl) : dE(di) indicates dwarf ellip- 
ticals with a certain, probable or possible underlying disk (i.e. 
showing spiral arms, edge-on disks and/or a bar) or other struc- 
tures (such as irregular central features (VCC21)); dE(bc) refers 
to galaxies with a blue center; dE to galaxies with no evident un- 
derlying structure. Four out of our 21 dwarf galaxies were not in 
the Lisker et al. sample (NGC3073, PGC 10072 17, PGC1 154903 
and VCC 1947), therefore we classified them as described in 
Section 6 attending only to their boxyness/diskyness. Column 6 
gives the Virgo substructure to which the galaxy belongs, taken 
fr om GOLDMine Data base dGavazzi et alJ l2003) and defined as 
in iGavazzi et a i1 (fl999l) . Columns 7 and 8 refer to the observa- 
tional campaign (see Table [1} and the exposure time for each 
setting, which allowed us to get a typical signal-to-noise raticQ 
for the central spectra of ~60 A -1 , enough to obtain reliable cen- 
tral kinematics. The galaxies were observed along their major 
axis. Their position angles (PA) from the HyperLEDA Database 
dPaturel et al.l l2003) are giv en in column 9. Colum n 10 gives the 
Galactic colour excess from lSchlegel et al.l ([1998). 

A to tal of 37 B to M stars in com mon with the MILES 
library dS anchez-Blazq uez et alJ 120061) and the CaT library 
dCenarro et alJ 120011) were observed to flux-calibrate our data 
and to use them as templates for velocity dispersion measure- 
ments. 



1 The S/N per A -1 is obtained dividing the S/N per pixel by the 
square root of the spatial scale along the slit. These measurement is 
therefore independent of the instrument used. 
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Run 1 


Run 2 


Run 3 


Date 


Dec.24-27 2005 


Jan.21-Zo ZUU / 


Feb. 10- 12 2007 


Tplpsmnp 


WHT 4.2m 


TNT 2 5m 


WHT 4.2 


Spectrograph 


ISIS 


IDS 


ISIS 


Detector 


EEV12(blue) 


Marconi(red) 


EEV10 


EEV12(blue) 


RedPlus(red) 


Grating 


R1200B(blue) 


R600R (red) 


R1200B 


R1200B(blue) 


R600R(red) 


Wavelength range 1 (A) 


3500-4300 


5500-6700 


3700-4790 


3500-4300 


5500-6700 


Wavelength range 2 (A) 


4100-4900 


7750-8950 


4600-5690 


4100-4900 


7750-8950 


Wavelength range 3 (A) 


4800-5600 










Dispersion (A/pixel) 


0.44 


0.87 


0.48 


0.44 


0.97 


Spectral Resolution (FWHM, A) 


1.56 


3.22 


1.80 


1.56 


3.23 


Instrumental Resolution (km s _1 ) 


40 


58 


46 


40 


58 


Spatial scale ("/pix) 


0.40 


0.44 


0.40 


0.40 


0.44 


Slit width (") 


1.95 


1.95 


1.94 


1.91 


1.91 



Table 2. The observed galaxies. 



Galaxy 


Other name 


RA(J2000) 


Dec.(J2000) 


Type 


Env. 


Run 


texp 


PA 


E(B-V) 






(h:m:s) 


(•:•:") 








(sec) 


C) 


(mag) 


M32 


NGC 221 


00:42:41.84 


+40:51:57.4 


cE2 


M31 group 


1 


1350 


170 


0.062 


NGC 3073 


UGC 05374 


10:00:52.10 


+55:37:08.0 


dE(di) 


Field 


1 


3200 


120 


0.010 


PGC 1007217 


2MASX J024135 14-08 10243 


02:41:35.8 


-08:10:24.8 


dE(di) 


Field 


1 


3600 


126 


0.024 


PGC 1154903 


2MASX J02420036+0000531 


02:42:00.3 


+00:00:52.3 


dE 


Field 


1 


3600 


126 


0.031 


VCC21 


IC 3025 


12:10:23.14 


+ 10:11:18.9 


dE(di,bc) 


Virgo N Cloud 


2 


3600 


99 


0.021 


VCC 308 


IC 3131 


12:18:50.77 


+07:51:41.3 


dE(di,bc) 


Virgo B Cluster 


3 


2400 


109 


0.021 


VCC 397 


CGCG 042-031 


12:20:12.25 


+06:37:23.6 


dE(di) 


Virgo B Cluster 


3 


3600 


133 


0.020 


VCC 523 


NGC 4306 


12:22:04.13 


+ 12:47:15.1 


dE(di) 


Virgo A Cluster 


1 


3400 


144 


0.044 


VCC 856 


IC 3328 


12:25:57.93 


+ 10:03:13.8 


dE(di) 


Virgo B Cluster 


2 


2740 


72 


0.024 


VCC 917 


IC 3344 


12:26:32.40 


+ 13:34:43.8 


dE 


Virgo A Cluster 


3 


3600 


57 


0.032 


VCC 990 


IC 3369 


12:27:16.91 


+ 16:01:28.4 


dE(di) 


Virgo A Cluster 


2 


3000 


135 


0.028 


VCC 1087 


IC 3381 


12:28:17.88 


+ 11:47:23.7 


dE 


Virgo A Cluster 


3 


3600 


106 


0.026 


VCC 1122 


IC 3393 


12:28:41.74 


+ 12:54:57.3 


dE 


Virgo A Cluster 


3 


3600 


132 


0.021 


VCC 1183 


IC 3413 


12:29:22.49 


+ 11:26:01.8 


dE(di) 


Virgo A Cluster 


2 


3600 


144 


0.031 


VCC 1261 


NGC 4482 


12:30:10.35 


+ 10:46:46.3 


dE 


Virgo A Cluster 


2 


6930 


133 


0.028 


VCC 1431 


IC 3470 


12:32:23.39 


+ 11:15:47.4 


dE 


Virgo A Cluster 


2 


3000 


135 


0.054 


VCC 1549 


IC 3510 


12:34:14.85 


+ 11:04:18.1 


dE 


Virgo A Cluster 


2 


3300 


13 


0.030 


VCC 1695 


IC 3586 


12:36:54.79 


+ 12:31:12.3 


dE(di) 


Virgo A Cluster 


3 


3600 


39 


0.045 


VCC 1861 


IC 3652 


12:40:58.60 


+ 11:11:04.1 


dE 


Virgo E Cloud 


3 


3600 


109 


0.030 


VCC 1910 


IC 809 


12:42:08.68 


+ 11:45:15.9 


dE(di) 


Virgo E Cloud 


2 


3800 


135 


0.030 


VCC 1912 


IC810 


12:42:09.12 


+ 12:35:48.8 


dE(bc) 


Virgo E Cloud 


2 


3600 


166 


0.032 


VCC 1947 


CGCG 043-003 


12:42:56.36 


+03:40:35.6 


dE(di) 


Virgo S Cloud 


2 


3060 


126 


0.027 



4. Data reduction 

The data reduction was performed with ReD^E dCardiell 19991) . 
a package specially designed to reduce long-slit spectroscopy 
with particular attention to the treatment of errors. This package 
is ideal for treating in parallel the data and error frames, pro- 
ducing an error spectrum associated with each individual data 
spectrum, which means that the errors are controlled at all times. 

Due to the similar instrumental configurations used on all ob- 
serving runs, the reduction process for both telescopes was the 
same. The standard procedure for long-slit spectroscopy data re- 
duction consists of bias and dark current subtraction, flat-fielding 
(using observations of tungsten lamps and twilight sky to cor- 
rect for high and low frequency variations respectively), cos- 
mic ray cleaning, C-distortion correction, wavelength calibra- 
tion, S-distortion correction, sky subtraction, atmospheric and 
interstellar extinction correction and flux calibration. We give 
below some comments on steps of particular importance: 



Flat-fielding. The flat-fielding correction is a delicate step 
at near infrared wavelengths due to the fringing effects. In the 
first run, the Marconi CCD suffered from significant fringing 
that varied with the telescope position. Since complete removal 
of the fringing in run 1 was not possible, we did not use the 
red Marconi-CCD data to determine the galaxy kinematics. The 
fringing produced by RedPlus, the new CCD optimised to avoid 
these patterns, was much lower, with an amplitude of only ~ 1 % 
independent of position. 

Wavelength calibration. The wavelength calibration was per- 
formed using between 65-100 arc lines depending on the instru- 
mental configuration. They were fitted with a 5' h order polyno- 
mial that led to a typical RMS dispersion of 0. 1-0.25 A. 

S-distortion, alignment of the spectra. During the spectro- 
scopic observations, the galaxies were not perfectly aligned with 
the rows of the detector. This effect is crucial when measuring 
gradients of any type (rotation curves, velocity dispersion pro- 
files or line-strength indices). The correction of this effect was 
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performed using a routine that found the position of the galaxy 
center as a function of wavelength, fitted all these positions with 
a low order polynomial and straightened the spectra using that 
polynomial. This alignment was done with a technique that min- 
imised the errors due to the discretization of the signal. This 
technique consists of adopting a more realistic distribution of 
the light in each pixel than just assuming it to be constant. To 
achieve this the signal in each pixel is fitted with a second or- 
der polynomial using the available information in the adjacent 
pixels. 

Sky subtraction. The sky subtraction is critical for studies 
where the spectra are analysed at light levels corresponding to 
only a few per cent of the sky signal, as in our case. For each 
galaxy observation a sky image was generated fitting the data 
at each wavelength with a first order polynomial in regions at 
both sides of the galaxy close to the ends of the slit (which has 
a length of 3.7 arcmin on the WHT and 3.3 arcmin on the INT). 
This was possible since for all targets except M32 the galaxy 
filled only a small region of the slit, so this synthetic sky im- 
age was free from contamination from the galaxy. For M32, we 
observed a separate sky frame moving the telescope from the 
coordinates of the galaxy to a position Aa =-416" (West), AS =- 
459" (South) far enough from M32 to avoid its light but with the 
same level of contamination from M3 1 . 

Extinction correction. Atmospheric extinc- 
tion was calculated using the extinction curve 
for El Roque de los Muchachos Observatory 
(www.ing.iac .es /Astronomy / observing/manuals /ps/ tech jiotes / 
tn031.pdf). T he Galactic extinc tion was corrected using 
the curve of jFitznatrick (1999) and the reddening from 
ISchlegel et ail dl998|) listed in Table [2] 

Flux calibration. The relative flux calibration of the spec- 
tra was perform ed using the observed stars in c ommon with the 
MILES library (ISanchez-Blazquez et al.l 120061) for the optical 
spectra, and with the CaT library ( Cenarro et al.l l200ll) for the 
near infrared. For each observed star we obtained a flux calibra- 
tion curve. All of them were averaged to obtain one unique flux 
curve for each run and instrumental configuration. The devia- 
tions of each flux calibration curve from the averaged one were 
introduced as uncertainties in the error spectra. The typical de- 
viation was of 2% reaching ~1% in the first and last ~150 A of 
each setup spectra where the noise is the highest. 



5. Measurement of the kinematic parameters 

The stellar kinematics of galaxies (radial velocities and velocity 
dispersions) were cal culated using the routine MOVEL included 
in REDmE package (Cardiel [l999l) . T his routine is based on the 
Fourier quotient method described bv lSargent & Turner! (Il977l) 
and refined with the OPTEMA algorithm dGonzalezll 1993b that 
allows us to overcome the typical template mismatch problem. 
In order to do this, a number of stars of different spectral types 
and luminosity classes were introduced in the program to cre- 
ate a model galaxy. These stars were of spectral type B9, AO, 
A3V, GO, G2III, G5III, G8III, G9III, KOIII, KOI, K2III, K3III, 
MOIII and M2III. The model galaxy was created and processed 
in parallel to the galaxy spectrum. To build the model galaxy all 
the template spectra were scaled, shifted and broadened accord- 
ing to a first guess of y (mean line-strength), v (radial velocity) 
and cr (velocity dispersion). Then the algorithm looked for the 
linear combination of these template stars that best matched the 
observed galaxy spectrum. The best linear combination of ob- 
served stars was chosen as the one that minimises the residuals 



between the galaxy spectrum and the broadened optimal tem- 
plate. This provided a first model galaxy with a first kinematic 
output (y, v and cr). This model galaxy was then improved us- 
ing this new guess of kinematic parameters. The process was 
iterated until it converged. The emission lines, found only for 
the field dwarf galaxies, and some large sky line residuals, only 
present in some cases, were masked, so that the program did not 
use them for the minimisation of the residuals. 

To minimize template mismatch effects, it is essential to use 
as templates a variety of spectral types and luminosity classes 
which are representative of the stellar population of the observed 
galaxy; as we will discuss in Section l5.3.11 small differences in 
ages and metallicities could lead to a partial fit of the strongest 
lines and therefore affect the derived velocity dispersion of the 
galaxy. 

It is also important to check whether the observed stars were 
filling the slit during the observation. If they were not, the instru- 
mental profile would not affect them in the same way as in the 
galaxies, and as a consequence, the cr that one would measure 

for the galaxy would be Jo~f nst + o~^ al , where cr ga i is the intrin- 
sic velocity dispersion of the galaxy and cr,„ Jf , in this case, is the 
quadratic instrumental difference between the galaxies and the 
stars. 

To correct for this effect the physical slit widtfQ is required. 
We calculated it from the spatial scale and the FWHM (in pixels) 
of the arc lines, which illuminate homogeneously the slit. To see 
whether the stars were filling the slit completely we checked that 
the FWHM of their spatial profile was larger than the physical 
slit width. If this was not the case, the spatial profile of the star 
was broadened accordingly. Although this introduced some un- 
certainties, the data quality improved by making this correction. 
The value of cr measured and adopted in this work is thus the 
intrinsic velocity dispersion of the galaxy corrected for possible 
instrumental effects. 

Figure [TJ shows a typical fit of the observed central spectral 
of a galaxy and the corresponding optimal template broadened 
with a gaussian with the derived dynamical parameters. The er- 
rors in velocity and cr were computed through Monte-Carlo sim- 
ulations, repeating the whole process (including the derivation of 
the optimal template) for 100 simulated galaxy spectra created 
using the error spectra obtained during the reduction process. 
The observed and simulated spectra perfectly match. 

5. 1 . Kinematic parameter profiles 

To measure kinematic gradients it is important to determine the 
minimum S /N needed to measure reliable radial velocities and 
velocity dispersions. To do that we have designed and carried out 
a test exercise based on Monte-Carlo simulations that constraints 
the errors and systematic effects in the measurement of radial 
velocities and velocity dispersions on fake galaxy spectra with 
known input kinematic parameters and different ages, metallici- 
ties and S/N ratios. This can be outlined in the following steps: 
i) from the simple stellar population models from PEGASE.HR 
(FWHM- 0.55 A: lLe Borgne et alJ (I2003I) ) we selected a subset 
of 9 model spectra (3 ages and 3 metallicities) representative of 
quiescent dwarf galaxies: ages 1 Gyr, 4 Gyr, 10 Gyr; metallic- 
ities of Z = 0.0, Z = -0.4, Z = -0.7. The spectral resolution 
of these models are needed since in our bluer configuration the 
resolution is ~ 1.6 A; ii) Taking into account our instrumental 



2 During the observation an indicative width of the slit is selected by 
the user through a web interface. 
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4200 



4400 4600 
Wavelength (A) 



4800 



Fig. 1. Example of the fit computed by MOVEL. Upper panel: in 
black, the central spectrum of VCC523. The spectrum has been 
continuum subtracted and normalised. The optimal template that 
fits the galaxy is shown in red, a linear combination with differ- 
ent weights of the stars observed with the same configuration 
as the galaxy. Bottom panel: residuals of the galaxy-composite 
template fit. 
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4100-4900A 



in 
d 



in 



~ o-=20km/s □ !j=!30km/s a' a=40km/d 

A cr=60km/s 



a=bukm/s . 




t=4Gyr Z=-0.4 
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Fig. 2. Simulations to study the minimum S/N ratio to obtain 
reliable measurements of the velocity dispersion. It is plotted 
Acr/cr (defined as (cr„ - cr,)/cr,), the relative error introduced 
when measuring the cr of a galaxy as a function of S/N ratio. 
Different colours show different cr for the simulated galaxies. 
Done assuming a stellar population of 4 Gyr and Z = -0.4. 



resolution, each model was broadened and redshifted to match 
a set of input velocity dispersions, cr,, (9 values between 20 km 
s and 60 km s in steps of 5 km s _I ) and radial velocities, 
v,, (800 km s _1 and 1500 km s _1 , typical values of Virgo cluster 
members). This amounts a total of 162 model spectra, iii) For 
each one of the above spectra we added different levels of ran- 
dom noise to match S/N ratios of 10, 15, 20, 25, 30 and 50, hence 
ending up with 972 model galaxy spectra, iv) For each simulated 
galaxy spectrum we run exactly the same MOVEL procedure 
as we did for our dE galaxy sample, using the same template 
stars and MOVEL parameters. 100 Monte-Carlo simulations for 
each model galaxy were carried out to get reliable errors of the 
derived kinematic parameters, cr and v„ obtained as the mean 
value of the 100 Monte-Carlo simulations in each case. Since 
the input kinematics cr,- and V; are set by construction, compar- 
isons and reliability analysis are immediate to perform. 

The above procedure was carried out for each instrumental 
configuration in each of the 5 different spectral regions. The re- 
sults obtained are shown in Figures [2] [3] and|4] 

After correcting for any systematic offsets in radial veloc- 
ity and velocity dispersion due to small intrinsic differences be- 
tween PEGASE.HR models and our observed stars, we have 
analysed the simulations looking at the relative differences be- 
tween the measured values and the parameters introduced in the 
simulated galaxies (Av/v = (v — Vj)/v; and Acr/cr = cr -cri)/cri). 
Figures |2] [3] and |4] show these differences as a function of S/N 
in the wavelength range 4100-4900 A, a range in common be- 
tween the WHT and INT observations and where lines as impor- 
tant as the G-band are located. The error bars in these 3 Figures 
show the relative uncertainties obtained by MOVEL as the RMS 
scatter resulting from the 100 Monte-Carlo simulations for each 
model galaxy. 

In Figure |2] we study the influence of the S/N ratio and the 
instrumental resolution on the measurement of the velocity dis- 
persion of a galaxy. Each point represents a galaxy of similar 
stellar populations (age 4 Gyr and metallicity -0.4) but differ- 
ent velocity dispersion (from 20 to 60 km s _1 ). As expected, the 
errors increase dramatically at the lowest S/N ratios. For low 



S/N ratios (S/N=10) offsets are found even for galaxies with 
high cr's, so the velocity dispersions derived at this S/N cannot 
be trusted. On the contrary, for S /N ratios higher than or equal 
to 15 we do not find statistically significant offsets, with the ex- 
ception of measurements below half the instrumental resolution 
(cr = 20 km s _1 ) where special care must be taken. Only for S /N 
larger than 20 the measured cr can be fully trusted for velocity 
dispersions as low as half the instrumental resolution. 

Figure [3]presents the influence of the stellar populations on 
the measurement of the velocity dispersion of a galaxy. In the 
upper panel we show the effect of the age on a dwarf galaxy of 
cr = 40 km s _1 and Z = -0.4, while in the lower panel we show 
the effect of the metallicity for a galaxy with cr =40 km s _1 and 
4 Gyr old. In this case, although the errors barely depend on 
metallicity, the age-dependence is crucial, and for populations 
as young as 1 Gyr the cr measurements are underestimated for a 
S/N ratio below 15. 

In Figure [4] we analyse the effect of the stellar populations 
on the computation of the radial velocity. Neither a change in 
the velocity dispersion nor in the metallicity have appreciable 
effects on this variable. Only the age of the stellar population 
have appreciable influence on the determination of the radial ve- 
locity whenever the age is young, around 1 Gyr. The uncertainty 
induced by age variations, however, is small, <1%, thus rotation 
velocities can be accurately measured down to S/N~10. 

In Figure|5]we plot the expected RMS error (i.e. RMS of the 
average of all the simulations) as a function of S /N ratio in order 
to have a statistical estimate of the uncertainty. The average of all 
the simulations (corresponding to different velocity dispersions, 
ages and metallicities) for each S/N, takes into account that for a 
target galaxy we do not have a priori information about either its 
velocity dispersion or the parameters of the stellar population. 

Figures [2]|5] clearly show that data with S/N ratios below 15 
might induce errors as large as 22% in the determination of cr, in 
particular for small velocity dispersions (cr ~20 km s~'), while 
only 0.4% for radial velocities. 

All these tests have been computed for the different wave- 
length ranges covered by our survey and the results obtained are 
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Fig. 3. Simulations to study the dependence of the stellar popula- 
tion of the galaxy on the measurement of its velocity dispersion. 
We plot the relative offsets found between the measured cr and 
the velocity dispersion introduced in the simulated galaxy as a 
function of the S/N ratio. In both panels a typical cr of 40 km 
s _1 has been considered. In the upper panel we fix the metallic - 
ity of the galaxy to -0.4 and study the influence of the age of the 
stellar population. In the lower panel the parameter fixed is the 
age to 4 Gyr, and we analyse the influence of the metallicity on 
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Fig. 4. Simulations to study the influence of the stellar popula- 
tions and the S/N ratio on the measurement of the radial veloci- 
ties. We plot the relative uncertainty when measuring the radial 
velocities (Av/v, defined as (v a - v,)/v,) as a function of the S/N 
ratio. Only the results for the radial velocity v =1500 km s _I 
are shown because the offsets and errors found are independent 
of the radial velocity of the galaxy. The results plotted are for a 
typical dwarf galaxy with velocity dispersion of 40 km s _1 and 
metallicity -0.4. The velocity dispersion does not have any ef- 
fect on the measurement of radial velocities nor the metallicities. 
Note the different scale in the y-axis compared to Figures [2] and 
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Fig. 5. Total scatter in the differences found for the 972 simula- 
tions computed as a function of the S/N ratio. In this Figure all 
models for all stellar populations and velocity dispersion param- 
eters have been used. 



rather similar. For the red arm of ISIS, where the instrumental 
resolution is larger, we obtain similar results as those shown in 
the blue arm (Figures[2]|5]) but with slightly larger uncertainties. 

These simulations show that radial velocities can be com- 
puted with spectra of S/N as low as 10 given that the uncertainty 
is always below 1%. However, in the study of velocity disper- 
sions, we must discard any measurement with S/N below 15 be- 
cause cr is, in this case, highly dependent on age and not reliable 
for cr as low as half the instrumental resolution. 

When running the MOVEL algorithm as a function of galaxy 
radius, we fitted the optimal template at every radius, rather than 
using the central optimal template, in order to improve the fit. 
As a result, the optimal templates turn out to be radial depen- 
dent. The differences, however, are not very large, because in 
the linear combination of templates, G-stars always contribute 
with the highest weight. 

Due to the different instrumental configurations used in the 
observation campaigns (see Table Q]), more than one kinematic 
profile per galaxy was obtained. These profiles, consistent within 
the errors, were averaged to produce a single, high S /N profile 
per galaxy (see Figure[6]). 

The recessional velocity of each single galaxy, removed for 
the determination of the rotation curve (Figure [6), has been de- 
termined by averaging, with a weighted mean, the recessional 
velocity measured in each single position along the radius. This 
improved technique for measuring recessional velocities (listed 
on Table|U) can be applied since the rotation curves are symmet- 
ric alQ 

Table [3] gives an example of the tables electronically avail- 
able with the values of the kinematic profiles. 

5.2. Central velocity dispersion and maximum rotational 
velocity 

To compute the central velocity dispersion (cr) we shifted all the 
spectra to the same wavelength scale using the rotation curves 
displayed in Figure [6] and we coadded all the individual spectra 



3 Note that galaxies as VCC856, those with the poorest quality, have 
a non-zero central velocity. 
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Fig. 6. Kinematic profiles of the galaxy sample. Each diagram shows in the left upper panel the folded rotation curve of the galaxy 
and in the left bottom panel the folded velocity dispersion profile. The different sides of the galaxy are indicated with red squares 
and black dots. On the upper x-axis the radius is given as a fraction of the effective radius (R e ff) of each galaxy in the I band (see 
Section |6). The purple open squares show the points used to calculate the maximum rotation for each galaxy and the purple line 
indicates this v max - The dashed line in the velocity dispersion profiles indicate the central cr computed up to the R e ff (see Table 0). 
In the right panels the not-folded kinematical profiles are plotted. 
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Table 3. Kinematic profiles for VCC990. Table 4. Kinematic parameters. 



Rv (") 


v (km s 1 ) 


R<r (") 


cr (km 


s- 1 ) 


-5.43 


1.8 ±5.5 


-4.98 


36.0 ± 


6.2 


-3.56 


16.8± 6.4 


-2.36 


47.2 ± 


6.0 


-2.58 


9.8± 6.3 


-1.38 


39.2 ± 


5.4 


-2.00 


10.5± 6.4 


-0.80 


35.8 ± 


4.6 


-1.60 


1.7+ 6.2 


-0.40 


32.5 ± 


3.5 


-1.20 


-6.3± 4.4 


0.00 


42.2 ± 


3.1 


-0.80 


3.9 ±2.7 


0.40 


41.3 ± 


2.8 


-0.40 


0.1 ± 2.3 


0.80 


39.7 ± 


3.6 


0.00 


6.8 ± 2.0 


1.38 


37.7 ± 


3.9 


0.40 


2.2 ± 2.4 


2.36 


36.1 ± 


8.5 


0.80 


1.0 ± 3.1 


4.97 


37.4 ± 


5.6 


1.20 


-5.2± 5.2 








1.60 


-11.8± 5.6 








2.00 


-25.3± 6.4 








2.58 


-17.1± 7.2 








3.56 


-29.6± 5.7 








5.42 


-25.3± 6.2 









NOTES: Column 1: radius for the rotation speed profile. Column 2: ro- 
tation velocities. Column 3: radius for the velocity dispersion profile. 
Column 4: velocity dispersions. All the kinematic profiles are electron- 
ically available. 



up to one effective radius. The typical S /N ratio for the spectrum 
where the central cr is computed is ~60 A . These results are 
shown in Table [4] 

The maximum rotational velocity (v max ) was calculated as 
the weighted average of the two highest velocities along the 
major axis on both sides of the galaxy at the same radius (for 
non-symmetric profiles at least three values were required). As 
a consequence values with larger errors weight less than those 
with smaller errors. For each galaxy these values are presented 
in Figure|6]as purple squares. We show in Appendix[A]that, given 
the uncertainty, all galaxies with v max < 9 km s _I can be consid- 
ered non-rotators. 

The ratio between these two kinematic measurements, the 
maximum rotation velocity v max and the velocity dispersion cr, 
is called anisotropy parameter, v max l<T, and it is used to study 
the rotational/pressure support of the galaxies. In Table [4] we 
show (Vmaxjcf)*, the anisotropy parameter corrected from the 
inclination. This correction is done following the expression 
(v m ax/o~Y - ^'™i(\- e ) ' wnere e is the ellipticity. Note that for 
those galaxies with ellipticity close to zero, no correction can 
be done because they are nearly face on. We choose a conserva- 
tive value of (Vmax/cr)* = 0.8 as the limit between pressure and 
rotationally supported systems in order to include those objects 
that, within the errors, are consistent with being flattened by ro- 
tation. This assumption is justified by the fact that the measured 
Vmax is a lower limit since the rotation curves are still rising. 

5.3. Comparison to the literature 

Displaying simultaneously the kinematic profiles measured in 
this work with those of other authors (Figure |7), one sees that 
the radial extent of the kinematic curves varies from one work to 
another. In addition, the offsets found in the velocity dispersions 
are not always consistent within the errors (Figure[8]). These two 
differences are important, because different radial extents lead to 
different maximum rotation velocities and offsets in the velocity 
dispersion profiles lead to different central cr values. 



Galaxy 


cr (km s 1 ) 


Vmax (km s ) 


(Vmax/v)* 


Vrad (km s ') 


PGC 10072 17 


35.2± 


0.9 


30.6± 5.0 


2.3 


± 0.4 


1592.3 


±3.3 


PGC 1154903 


23. 1± 


4.1 


8.6± 2.6 


0.5 


± 0.2 


1156.3 


±6.4 


NGC3073 


39.8± 


0.3 


17.1± 6.7 


1.1 


± 0.4 


1168.9 


±2.5 


VCC21 


26. 1± 


4.0 


18.4± 5.9 


0.9 


±0.3 


463.6 


t9.2 


VCC308 


31.7± 


1.2 


30.4± 8.6 


1.0 


±0.3 


1515.3 


± 3.2 


VCC397 


29.9± 


1.1 


52.0± 11.6 


2.5 


± 0.6 


2434.9 


±2.4 


VCC523 


45. 8± 


0.7 


39.6± 5.7 


1.5 


± 0.2 


1515.7 


± 2.6 


VCC856 


29.6± 


2.5 


9.7± 1.9 


1.0 


± 0.2 


980.2 


t 1.9 


VCC917 


31.4± 


1.4 


21.6+ 7.5 


0.8 


± 0.3 


1236.2 


± 2.4 


VCC990 


40.6± 


1.0 


26.3± 1.6 


0.9 


±0.1 


1691.5 


± 2.1 


VCC1087 


48.3± 


0.7 


7.1± 6.4 


0.2 


± 0.2 


644.6 


t3.2 


VCC1122 


37.2± 


0.8 


17.3± 7.7 


0.5 


± 0.2 


447.8 


t 2.7 


VCC1 183 


41. 3± 


1.2 


10. 1± Z.O 


0.5 


± 0.1 


1290.7 


±2.3 


VCC1261 


51.8± 


0.9 


13.9+ 5.2 


0.4 


±0.1 


1806.4 


± 2.1 


VCC1431 


54. 1± 


1.2 


7.0± 3.6 


0.1 


±0.1 


1472.8 


± 3.1 


VCC1549 


38.9± 


1.9 


5.2± 2.2 


0.3 


± 0.1 


1377.3 


±2.9 


VCC1695 


28.7± 


1.1 


14.3+ 3.1 


0.9 


± 0.2 


1706.0 


± 3.0 


VCC1861 


40.4± 


0.9 


6.3± 4.3 


0.2 


±0.1 


617.2 


t2.6 


VCC1910 


39.0± 


1.1 


6.0± 2.3 


0.4 


± 0.2 


178.6 


t2.2 


VCC1912 


37.1± 


1.0 


20.8± 4.5 


0.5 


± 0.1 


-110.5 


± 2.0 


VCC1947 


45.3± 


1.1 


28.3± 2.1 


1.1 


±0.1 


956.3 


t 1.8 



NOTES: Column 1: galaxy name. Column 2: central velocity disper- 
sions computed within the R e ff- Column 3: maximum rotation veloc- 
ities. The v^x adopted for VCC0856, V CC0990 and VCC11 83 have 
been measured in the rotation curves of IChilingarfanl j2009h due to 
their larger extent (see Section 5.3). Column 4: anisotropy parameter 
corrected from inclination. Column 5: mean radial velocity observed. 
Values in agreement with those of NED database. 



5.3.1. Velocity dispersions 

There are two factors that critically affect the measurements of 
the velocity dispersion of the galaxies: the not identical instru- 
mental cr for the stars and the galaxies and the spectral types of 
the stars used to fit the w idth of the galaxy lines. 

ICaldwell et al.l (I2003I) made their observations with a cr inst = 
100 km s _1 , too high to accurately measure velocity dispersions 
of typical dwarf galaxies, below 50 km s . 

Concerning the equality of the instrumen tal re s olution 
in the templates and in the galaxies IChilingarianl £2009) 
used simple stellar pop ulation (SSP) models of PEGASE .HR 
dLe Borgne et alJ |2003), based on ELODIE (Prugniel et al. 
120071) (R = 10000). In this case, the broadening of the SSP mod- 
els to cr ins , of the galaxies must be cautiously done because cr msl 
of the models is based on the mean value of the spectral resolu- 
tion of the library they are based on, therefore small differences 
after the broadening between the SSP models and the galaxies 
can arise. 

Co ntinu ing with the same effec t, in the works of Pedraz et al.l 
( 2002J) and Ivan Zee et al. (2004b) the stars observed were not 
filling the slit because they were not defocused and the seeing 
was sm alle r than their slit width . But in the case of iGeha et al.l 
d2003l) and iBeaslev et ail d2009l) . where their slit widths were 
0.75" and 1.0" respectively, it is possible that some stars were 
filling the slit thanks to the seeing. However, in none of these 
works a correction was made to assure that the instrumental pro- 
file in the stars was the same as in the galaxies. 

In reference to the templates used to perform cr, some of 
the authors mentioned above used only one star to fit the galaxy 
spectrum, and in such a case the fact that the template is not 
representative of the stellar population of the galaxy might lead 



E. Toloba et al.: Formation and evolution of dwarf early-type galaxies in the Virgo cluster 



9 



Table 5. Uncertainties introduced when different templates are 
used to calculate the velocity dispersion of a dwarf galaxy with 
different stellar populations. Acr/cr, defined in Figure [2] uses as 
cri 40kms~ 1 . 



I ^xlOO 

1 IT 


Age (Gyr) 


Z=+0.0 


Z=-0.4 


Z=-0.7 


Linear combination of B to M stars 


1 


15% 


19% 


1% 


4 


8% 


18% 


19% 


10 


7% 


14% 


23% 


G8III 


1 


32% 


54% 


71% 


4 


15% 


18% 


20% 


10 


17% 


13% 


11% 


Kim 


l 


29% 


47% 


59% 


4 


19% 


18% 


20% 


10 


17% 


13% 


11% 



to large errors. We have computed Monte-Carlo simulations to 
see the differences between fitting the galaxy spectrum with only 
one star and a linear combination of stars of spectral types from 
B to M, with different luminosity classes. The simulations con- 
sisted of a selection of PEGASE.HR models of 3 different ages 
(1 Gyr, 4 Gyr and 10 Gyr) and 3 different m etallicities (Z = 
0.0, Z = -0.4, Z = -0.7) with Salpeter IMF (iLe Borgne et al l 
120031) : a total of 9 m odels. The stars used as templ ates were 
from MILES library dSanchez-Blazquez et ail 120061) . First of 
all we checked that after broadening the models (originally at 
FWHM=0.55 A) to the MILES resolution (FWHM=2.3 A) we 
obtained cr — km s _1 when running MOVEL, thereby showing 
that there was no zero point offset. Secondly, we broadened the 
models to 40 km s to simulate dwarf galaxies of different stel- 
lar populations. And finally, we ran MOVE L using 3 d i fferen t 
kinds of templates: only one K1III star (a s inlGehaet all (2003) 
), only one G8III star (as in Ivan Zee et al.l (l2004bl) ). and a linear 
combination of B to M stars with different luminosity classes (as 
in this work). To measure cr we have masked the Balmer lines, 
especially those bluer than Hg to avoid possible problems due to 
emission lines. Our simulations show that a linear combination 
of different stars is the most accurate method to obtain the veloc- 
ity dispersion of the galaxies, never finding an error above 25%, 
independent of the stellar population considered. Note that if a 
young population dominates the light of the galaxy (for ages of 
1 Gyr and below), and a single G or K star is used as a template, 
an error up to 70% can be done. When a single star is used as 
template, a dependence on metallicity for young populations (1 
or 4 Gyr) is also found, in the sense that decreasing the metal- 
licity increases the uncertainty. This dependence is likely due to 
offsets introduced by the method employed to compute cr. The 
results obtained have been summarised in Table [5] 

5.3.2. Rotation curves 

Different criteria have been used in the literature to measure 
the maximum rotation velocity. The main difficulty here is to 
have an extended rotation curve that reaches a clear plateau 
where the maximum rotation can be measured. As this is not 
so easy for dEs, an objective criteria, independent of the shape 
of the rotation curve in each case, must be adopted. The differ- 
ent criteria used by the various authors have led to maximum 
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Fig. 8. Comparison between the velocity dispersions measured 
in this work versus those of other authors. The colours and sym- 
bols are the same as in Figure [7] We add a comparison with 
ICaldwell et al.l d2003l) (solid red points), who only measured 
central values and no kinematical profiles. 



rotational velocities that are nevertheless in the majority of the 
cases n early consistent with in the errors (Figure|9]l- The criterion 
used by Pedra z et alj (|2002|) is the same as the one adopted here. 
Although we achie ved a radial extent of 23" for VCC1 122 and 
Pedr az et"aT1 J2002) only 8", the latter value was enough to reach 
the flat part of the rotation curve, and, as a consequence, both 
measurements of the maximum rotatio n are identical. The expla- 
nation for the differences found with Simien &PrugnieH (2002) 
are mainly based on the fact that only two points were considered 
to ob tain the maximum rotation in that paper. The differences 
with iGeha et all (120031) are due to the limitation of their data 
to the core of the galaxies, never re aching radii larger tha n 6" 
(see F igure |7V The differences with Ivan Zee et al.l d2004bl) and 
Chilin garian (2009) are related t o a different radial extent of the 
rotation curves. Note that Chilingarian (2009) does not calculate 
the maximum rotation, but we have applied our c riterion to his 
rotatio n curves. Finally, the differences found with Beasley et al. 
(2009) are due to the fact that v max is obtained from the anal- 
ysis of Globular Clusters located up to ~ !R e ff- When their 
kinematic determined from the stellar component using long slit 
spectroscopy along the major axis of the galaxy is compared to 
our data, the agreement is evident (Figure |7). The maximum 
rotation values adopted for the analysis (see Table [U are our 
own values, except for VCC856, VCC9 90 and VCC1 183, where 
the data come from Chilingarian ( 2009 ) since he obtained la rger 
radii than us. Note that the values from lBeaslev et al] (2009) can 
not be adopted here because our work is dedicated to the analysis 

of the stellar component of dEs. 

The comparison with iBeaslev et alJ ([2009), who finds rota- 
tion speeds much larger at !R e ff (100 and 50 km s _1 higher for 
VCC1087 and VCC 1261, respectively) than our data at the R eff , 
indicate that the rotation curves of these galaxies are still rising. 



6. Photometric parameters 

In order to make a complete analysis of the kinematics, compar- 
ison with some photometric parameters is needed. For our study, 
we require I-band (Johnson-Cousins) total magnitudes and opti- 
cal rad ii (R ODt , radius contai ning 83% of the total I-band lumi- 
nosity (Catinella et al. 2006)) to study the shape of the rotation 



10 



E. Toloba et al.: Formation and evolution of dwarf early-type galaxies in the Virgo cluster 



V) 



6 

In 



O 
00 



o 

CO 



o 



o 

CN 



o 

E o 



Pedraz et 'al. 2002 
Simien Sc Prugniel 2002 
Geha et al. 2003 
Van Zee et al. 2004 
Chilingarian 2009 
Beasley et al. 2009 




to 







20 40 
(km/s) This Work 



Fig. 9. Comparison between the maximum rotational velocity 
measured in this work and those measured by other authors. 
Sym bols and colours are the same as in Figure [7] VCC1261 
from Beasley et al. (2009) is nearly outside the plot due to its 
enormous rotation: 105 +44 km s -1 . 



curves. Effective radii (R e ff, radius containing 50% of the to- 
tal light) is needed to measure the extent of the radial profiles 
in physical units of the galaxy. Ellipticities (e) are needed to 
make the appropriate corrections due to inclination. A param- 
eter to measure the boxyness/diskyness of the isophotes (C4) is 
also required to study the possible late-type origin of these dwarf 
early-type galaxies. 

All these parameters hav e been obta i ned fr om ;-band Sloan 
Digita l Sky Survey (SPSS. lYork et al.l (I2OOOI) ') data release 6 
(DR6, Adelman-Mc Carthv et al.l d2008b ') photometry. They have 
been calculated using the IRArfj task ellipse as described in 
Appendix iBl The transformation from /-band (SDSS) to /-band 
(Joh nson-Cousins) has bee n done assuming mi = m, -0.52+0.01 
mag dFukugita etal .1119951) . 

6.1. C4: boxyness/diskyness parameter 

The boxyness/diskyness (C4) parameter measures the deviations 
of the isophotes from a perfect ellipse. If C4 > the isophotes 
are disky, indicating that some disk substructure is present in the 
galaxy, and when C4 < the isophotes are boxy dCarterl [l978: 
Kormendy & B enderll996l) . This parameter is independent from 
the surface brightness profile of the galaxy. The C4 parameter, 
determined for our galaxies as described in Appendix|Dl is pro- 
vided in Table [6] 

Our C4 classification for disky isopho tes agrees in gene ral 
with the morphological classification of Lisker et al. (2006a|) as 
can be seen in Figure [TOl Note that lLiskeretalJ (l2006a) clas- 
sified a galaxy as disky when disk features (spiral arms, edge- 
on disks, or bars) were detected after subtracting an axisym- 
metric light distribution from the original image or after un- 
sharp masking. In this Figure we can see that the red dots, 
th ose galaxies classifie d as being without underlying structures 
in lLisker et al.l (l2006al) . are grouped around C4 < (boxy or el- 
liptical isophotes), while the blue dots (galaxies with some disky 



4 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under the cooperative agreement with 
the National Science Foundation. 
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Fig. 10. Correlation between the a nisotropic p a rameter 
(Vmax/cr)* and C4XIOO. Colours indicate iLisker et"aU (|2()06b) 
classification. The galaxies with v max measured inside the cen- 
tral 6" are considered lower limits and are indicated with arrows. 
The horizontal dashed line at (v max /<ry =0.8 indicates a rough 
boundary between pressure and rotationally suppor ted galaxies. 
The blu e points are classified as having a disk in ILisker et"aTI 
(2006b), whereas for the red points no disk is found in that pa- 
per. See Appendix ID1 for a discussion. 



structure in lLisker et al.l d2006al) ) are all consistent with C4 > 
(disky isophotes). This justifies the use of the C4 parameter to 
detect the presence of an underlying disk. Three exceptions are 
found and discussed in Appendix ID1 

In this respect, it is important to underline the correlation 
between C4 an d the anisotropic par ameter evident in Figure [TO] 
first found bv iBender etal] (fl988) for more massive elliptical 
galaxies. 



7. Analysis 

The analysis presented in this work is primarily focused on the 
rotationally supported systems. Although the majority of the 
dwarf galaxies (15 out of 21) show some rotation (y max > 9 km 
s _1 , Table [H), only 11 a re rotationally supported ((v max /o-y > 
0.8) (Toloba et al. 2009). Here we try to understand whether the 
observed kinematic properties of the rotationally supported sys- 
tems are consistent with those of star forming galaxies of similar 
luminosity. 



7.1. Shape of the rotation curves 

ICatinellaetail(l2006l) made a systematic study of the shape of 
the rotation curves of late-type spiral galaxies a s a function 
of lum inosity based on the method described in iPersic etaD 
(1996]). They fitted the rotation c urves following the Polyex 
model dGiovanelli & Havnesll2002l) which has the form: 



V/^(r) = Vo(l-^' /rre )(l + ^) 



(1) 



This analytical function depends on 3 parameters: V > fpE an d a > 
which represent the amplitude, the exponential scale of the inner 
region and the slope of the outer part of the rot ation curve, re- 
spectively. The mean fitted rotation curves from ICatinella et all 
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Galaxy 


d (Mpc) 


M, (mag) 


6 


Ref 


r (") 


C 4 xl00 


Disk classification (L06b) 


PGC 10072 17 


20.23 


± 1.40 


-17.39 ± 


0.15 


0.13 


±0.05 


11.11 


± 0.02 


4.8 ( *» 


± 0.6 




— 


PGC 1154903 


15.00 


± 1.10 


-15.91 ± 


0.16 


0.34 


±0.02 


8.36: 


t0.04 


0.3 = 


t4.4 




— 


NGC3073 


33.73 : 


t 14.44 


-20.43 ± 


0.93 


0.14 


±0.02 


14.78 


±0.03 


3.2" 


± 0.2 




— 


VCC21 


16.74 


± 0.15 


-17.78 ± 


0.03 


0.36 


±0.03 


15.48 


± 0.11 


3JW 


±0.6 




Disk 


VCC308 


16.41 


± 0.32 


-18.78 ± 


0.05 


0.04 


±0.03 


19.22 


± 0.04 


0.2 = 


t2.0 




Disk 


VCC397 


16.41 


± 0.32 


-17.62 ± 


0.05 


0.33 


±0.03 


13.75 


±0.02 


5.1W 


±0.3 




Disk 


VCC523 


16.74 


± 0.15 


-19.16 ± 


0.03 


0.25 


±0.01 


20.90 


±0.03 


2.6<*> 


±0.3 




Disk 


VCC856 


16.83 


±0.46 


-18.49 ± 


0.06 


0.09 


±0.03 


15.82 


± 0.06 


-0.2 


± 1.7 




Disk 


VCC917 


16.74 


± 0.15 


-17.39 ± 


0.03 


0.41 


±0.02 


9.68: 


t0.04 


4.3 ( *> 


±0.2 




No Disk 


VCC990 


16.74 


± 0.15 


-18.27 ± 


0.03 


0.34 


±0.02 


9.73: 


t0.02 


4.1« 


±0.3 




Disk 


VCC1087 


16.67 


±0.46 


-18.94 ± 


0.06 


0.28 


±0.03 


22.74 


±0.05 


-0.9 


± 1.1 




No Disk 


VCC1122 


16.74 


± 0.15 


-17.94 ± 


0.03 


0.50 


±0.04 


14.06 


±0.05 


3.8 W 


± 0.2 




No Disk 


VCC1183 


16.74 


± 0.15 


-18.59 ± 


0.03 


0.22 


±0.12 


18.23 


±0.01 


1.7W 


±0.0 




Disk 


VCC1261 


18.11 


±0.50 


-19.39 ± 


0.06 


0.37 


±0.05 


22.07 


± 0.03 


1.5 W 


± 0.6 




No Disk 


VCC1431 


16.14 


±0.45 


-18.47 ± 


0.06 


0.03 


±0.01 


10.32 


±0.01 


0.2 = 


t 0.1 




No Disk 


VCC1549 


16.74 


± 0.15 


-18.18 ± 


0.03 


0.16 


±0.01 


13.56 


± 0.05 


-1.3 


±3.5 




No Disk 


VCC1695 


16.52 


±0.61 


-18.13 ± 


0.08 


0.22 


±0.05 


19.78 


±0.10 


-0.9 


±3.7 




Disk 


VCC1861 


16.14 


±0.45 


-18.57 ± 


0.06 


0.04 


±0.02 


18.52 


±0.04 


0.8 = 


t2.5 




No Disk 


VCC1910 


16.07 


±0.44 


-18.63 ± 


0.06 


0.14 


±0.04 


13.70 


±0.02 


-0.4 


± 1.0 




Disk 


VCC1912 


16.74 


± 0.15 


-18.62 ± 


0.03 


0.54 


±0.06 


23.34 


±0.02 


3.6W 


± 0.2 




No Disk 


VCC1947 


16.74 


± 0.15 


-18.46 ± 


0.03 


0.23 


±0.01 


10.70 


± 0.02 


4.2<*» 


±0.5 







NOT ES: Colum n 1: galaxy name. Column 2: distances in Mpc from Surface Br i ghtne ss Fluctuations (SBF) for individual Virgo galaxies from 
iMei et all ((2007) when available, or the mean A/B Cluster distance from Mei et al. ( 2007) for the rest of them (note that E, N and S Clouds are East, 
North and South areas of Cluster A). The distance for NGC3073 from SBF bv lTonrv et al.1 d2001l) and PGC1007217 and PGC1 154903 distances 
are from NED/IPAC Database derived from redshift with Hq =73 ± 5 km s _I Mpc~' (these two distances must be used cautiously). Column 3: 
Absolute magnitudes in /-band (Johnson-Cousins in AB system) converted from /-band measured in SDSS images using m ; = m, - 0.52 mag. 
Column 4: ellipticities from /-band (SDSS) images. The quoted errors indicate the RMS scatter in the ellipticity between 3" and the R e ff- Column 
5: Effective radius from i-band (SDSS) images. Column 6: Diskyness/Boxyness parameter from i-band (SDSS) images. The asterisks indicate C4 
measured as the maximum in the region 3"- 3R e ff if prominent disky features are found; in the rest of the cases the quoted values are the average 
in this same radial range and the errors the RMS scatter (see Appendix iDt. Column 7: Disk/No Disk classification by Lisker et al. (2006b) (L06b) 
(no values indicate that these galaxies are not included in their analysis). 



(2006) are normalised to the optical radius (R opl , radius contain- 
ing 83% of the total /-band luminosity), and the velocities are 
corrected from inclination. To compare them with our rotation- 
ally supported galaxies, we have calculated the inclinations as in 
iGiovanelli et all dl997ah : 



2 . d-e) 2 -^ 

COS I = 



1-95 



(2) 



where i is the inclination, e is the ellipticity and qo is a constant 
value that depends on the thickness of the disk. Here we assume 
qo = 0.3, a conservative value for dwarf galaxies shaped as thick 
disks dLisker et al.ll2007l) PI 

The comp a rison between the mean rotation curves of 
ICatinella et al.l (120061) and those of our rotationally supported 
objects must be done in the same luminosity regime since the 
derived parameters of the Polyex model are luminosity depen- 
dent. As the dwarf galaxies analysed in this work have magni- 
tudes below the minimum magnitude in ICatinella et al.l (2006) 
{Mi = -19.4), the reference Polyex model of low luminosity star 
forming systems has been determined extrapolating linearly the 
three faintest values of the Polyex parameters to M/ = -18.49, 
the mean /-band magnitude of the galaxies here analysed. The 
parameters used to construct this curve are Vo = 74.58 km s _1 , 
rpE = 0.35" and a = 0.03. 



5 For early-type spirals go = 0.2 ( IGiovanelli et alJ[l997al) . The most 
recent mea surements of the thickness o f dwarf galaxies gives q = 
0.3 - 0.35 ( Sanchez- Janss en et a D l2010h . The difference in v max after 
the correction for inclination between using q = 0.2 and q = 0.35 is 
of 4.4%, insignificant. 



In Figure [TTj we c ompare the mean fitted rotation curves of 
ICatinella et al. ( 20061) for late-type spirals (obtained from emis- 
sion lines, black solid curves) with the rotation curves of our ro- 
tationally supported dEs determined from absorption lines (grey 
symbols). Of the 1 1 rotationally supported dEs, only 7 have been 
considered for this analysis because 3 of them, VCC21, VCC917 
and NGC3073, have poor quality rotation curves (-p 24 larger 
than 25%), and VCC308 has e lower than 0.1, implying that the 
galaxy is nearly face on. The blue dots in Figure Q~T] show the 
median rotation curve of our dEs in bins of =0.1. The grey 
area contains the Icr deviation from this median value (68% of 
the values are inside this area). The blue dashed line is the ex- 
trapolated Polyex model for M/ = -18.49. 

It is evident from Figure QT| that our rotationally supported 
galaxies are characterised by rotation curves that are similar to 
those of late-type spiral galaxies of equal luminosity. Or, in other 
words, galaxies with similar rotation curves have similar abso- 
lute magnitudes despite their morphological type. We see that 
rotationally supported dEs dynamically behave like small late- 
type spiral galaxies. 



It is interesting to see that, despite their similar exponential 
radial light distribution, dEs have two different kinematic be- 
haviour, they can be either pressure or rotationally supported. 
Furthermore, rotationally supported dEs have rotation curves 
similar to those late-type spirals despite the fact that these lat- 
ter objects are gas dominated systems. 
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Fig. 11. The observed rotation curves of rotationally supported 
dEs (grey symbols) are compared to the mean rotation curves 
of lat e-type spiral galaxies (black solid and blue dashed lines) 
from ICatinella et aD d2006l) . Blue filled dots represent the me- 
dian observed rotation curve of rotationally supported dE in bins 
of r/R opt = 0.1. The last bin contains all data for r/R op , > 0.5. 
The grey area indicates rotation velocities within lcr from the 
median. 



7.2. Tully-Fisher relation 

Given the similarity in the kinematic properties of rotation sup- 
ported dEs with those of late-type spirals we expect that these 
systems follow the Tu lly-Fisher relation, as firstly proposed by 
Ivan Zee et al.l (l2004bl) . The Tully-Fisher relation is a typical 
scaling relation valid for star forming, rotating systems, link- 
ing the total luminosity to the maximal rotation velocity of the 
galaxy. 

In Figure fT2l we compare the Tully -Fisher relation for our 
(dark bluefj and Ivan Zee et alj d2004bl) (light blue) dEs to that 
of normal late -type galaxies of IGiovanelli et all (Il997bl) (grey 
symbols) and De Riicke et al.l (120071) (red dashed area), respec- 
tively. Figure [12] clearly shows that these rotationally supported 
dEs follow the Tull y-Fisher relation wit h a similar scatter as the 
normal spirals of D e Riicke et al.l (120071) and thus kinematically 
behave as late-type spirals. The v max of dEs plotted in Figure [T2l 
is probably a lower limit since it is generally measured where 
the rotation curve is still rising, as suggested by the kinematics 
of the globular clusters. It is thus conceivable that the agreement 
between the Tully-Fisher relation of rotationally supported dEs 
and late-type spirals of similar luminosity is even better than that 
depicted in Figure [T2l 
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Fig. 12. Tully-Fisher relation for 7 of our rotationally sup- 
ported dEs (in dark blue), the dEs from Ivan Zee et all d2004albl 
VZ04, in light blue), the n ormal spirals from [G iovanell i et all 
d!997bl G97, in grey) and iDe Riicke et al] d2007l DR07. red 
limited area). Absolu te magnitu d es of dEs have been obtained 
using distances from iMei et all {2007) (criterion described in 
Table |6}. For Giovanelli et al. data we use Hq = 73 km 
s 1 Mpc- 1 dMei et all d200l ). The arrows indicate lower lim- 
its of Vmax (those obtained in the inner 6"). Fits of the Tully- 
Fisher relation are indic ated in black f or the normal sp i rals of 
IGiovanelli etal] d!997bl) and in red for IDe Riicke etal] d2007l 
DR07). DR07 fit has been trans formed to I band using the 
colour-morphology relation from Fukugita et alj d 1995b and us- 
ing H = 73 km s' 1 Mpc" 1 and M /Q = 4.08 and M Bo = 5.48 
from lBinnev & Merrifieldl dl998h . 



where M press is the mass inferred from the velocity disper- 
sion after the contribution from rotation has been removed and 
M rot is the mass deduced by the intrinsic rotation velocity of 
the galaxies. M press in side the half-light radius is defined as in 
Cappellari etall d2006l) : 



M press 2.5G~ 1 0- 2 /?, 



580 



R eff 

km 2 s~ 2 /\ pc 



M, 



(4) 



The rotation curves of rotationally supported systems are 
characterised by an approximately constant gradient suggesting 
solid body rotation up to the R e ff- In this case M rot is given by 
the relation: 



7.3. Dark matter content 

The shape of the rotation curves gives information about the 
dark matter content and distribution of late-type galaxies (e.g., 
ICatinella et al.ll20Q6h . Similarly, <x can be used to measure the 
dark matter content of pressure supported systems. Following 
Beasl ev et al.l d2009l) we estimate the total dynamical mass of 
our sample galaxies using the relation: 



M press M ro i 



(3) 



6 No asymmetric d rift is applied neither to our dEs nor to those of 
van Zee et al. ( 2004b). This correction would increase v„ mx by 2.5 ± 0.9 
in those rotationally supported systems plotted in Figure |T2] 



R effV, 



1 

mux 



G 

R eff 

pc /\kmV 2 /\4.3 x 10 



1 



M 



(5) 



Dynamical mass-to-light ratios (T/, Table [7]) are then mea- 
sured using the /-band luminosities and equation [4] for pressure 
supported systems and the sum of equations|4]and|5]for rotation- 
ally supported objects Q Stellar mass-to-light ratios (TJ, Table 
13 are computed using the models of single stellar populations 



7 Note that this method to obtain M tc 
asymmetric drift. 



is equivalent to introducing the 
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Table 7. Dynamical and stellar mass-to-light ratios in I-band in 
solar units. 



Galaxy 


00) o 


(rp Q 


PGC 10072 17 


4.0 + 0.8 


1.9 


PGC 1154903 


1.8 ±3.1 


1.2 


NGC3073 


0.6 ± 0.4 


0.2 


VCC21 


1.5 + 0.4 


0.5 


VCC308 


1.5 + 0.3 


1.1 


VCC397 


4.8 ± 1.2 


0.9 


VCC523 


2.2 ± 0.2 


1.3 


VCC856 


1.1 ± 0.2 


2.0 


VCC917 


2.5 ± 0.4 


2.4 


VCC990 


1.9 + 0.1 


2.8 


VCC1087 


2.3 + 0.1 


2.4 


VCC1122 


2.4 + 0.1 


2.2 


VCC1183 


1.9 ± 0.1 


1.4 


VCC1261 


2.3 + 0.1 


1.4 


VCC1431 


2.3 + 0.1 


3.4 


VCC1549 


2.1 ± 0.2 


3.1 


VCC1695 


1.7 + 0.3 


1.2 


VCC1861 


1.6 + 0.1 


2.8 


VCC1910 


1.6 + 0.1 


3.2 


VCC1912 


2.1 ± 0.1 


0.6 


VCC1947 


1.8 + 0.1 


1.7 



NOTES : (T*) n is only indicati ve, the large errors in the stellar popu- 
lations (Michielsen et al. 2008) make the uncertainties of (Tp G of the 
same order as the values. 



(SSP ) oflVazdekiset al.1 d2010h and the ages and metalhcities 
from|4ichielsenetaL| d2008j). 

Figure [13]presents the relation between the dynamical mass- 
to-light ratio and the absolute I-band magnitude for our galax- 
ies (red and blue dots), the sam ple of classical elliptical galax- 
ies fr om Cappellari et al.l (120061 T j eans in I-band), the dEs from 
iGeha et al.l(l2002l) andl he Milky Way dwar f spheroidals (dSphs) 
fro mlWolf et alJd201(#l Consistently with lZaritskv et al.ld2006h 
and IWolf etal.1 d2010). we observe that the total mass-to-light 
ratio within the R e ff of dEs is the lower limit of the decreas- 
ing and in creasing T/ vs. lumino sity relations observed in gian t 
ellipticals dCappellari et al.ll2006l) and dSphs dWolf et al.l l2010). 
respectively. Dwarf early-type galaxies have on average T/ = 
2.00 ± 0.04 Ti , with a slightly higher dispersion in rotating 
systems (RMS = 0.06 Ti Q ) than in pressure supported systems 
(RMS = 0.04 1T I0 ) . 

Our sample of dEs have on average Md y „/M* = 1.6 ± 1.2 
(obtained as the ratio between the dynamical and stellar mass- 
to-light ratios, T/ and T*, respectively, Table 13, thus they are 
not domin ated by dark matte r within the R e f f (as p reviously sug- 
gested bv lGeha etaTI (120021) : iForbes et alT(l2008[)). consistent! 



with what is found is massive elliptical s (Cappella ri et al.ll200 



•6) 



witn wnat is found is massive ellipticals 
but contrary to dSphs dWolf et al.ll20 1 Ot) . 



8. Discussion 

How does this observational evidence compare with the different 
scenarios of galaxy formation? In a more general context we re- 
call that in the most recent hierarchical models of galaxy forma- 
tion only the most massi ve ellipticals have be en formed through 
major merging events (iDe Lucia et al.ll2006l) . The strong mor- 
phological segregation observed in high density environments 
(Sandageetai. 19851: iFerguson & Binggeiilll994l: iBlanton et al.1 



8 The T ; values of Wolf et al. (2010) have been converted to equation 
|4]for consistency. 
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Fig. 13. Dynamical mass-to-light ratio as a function of the ab- 
solute magnitude in I-band. Red and blue dots are our pres- 
sure and rotationally support ed dEs respectivel y. For compari- 
son open circles are dEs by Geh a et aT| d2002l G02), grey as- 

' 120 id V 



terisks are dSphs from IWolf et akl d2010i W10) and dark and 



light grey triangl es are slow and fast rotators respectively from 
Cappellari et al] d2006l C06). The transformation to I band 
has been performed using the colour-morphology relation by 
Fukugita et al.l d 19951) for Es. For the d warf galaxies we hav e 
used V - 1 = 1 .03 + 0.04 as calculated bv lvan Zee et alj d2004j) . 
Note that NGC3073 in not a dE given its high luminosity (M/ = 
-20.4 mag). 



120051) indicates that the cluster environment plays a major role 
in the formation of dEs. 

In lToloba et al.l d2009l) we have shown that rotationally sup- 
ported dEs, characterised by a disky structure, are preferentially 
located in the field and in the periphery of the cluster, while pres- 
sure supported systems are closer to the center. We also found 
that rotationally supported dEs have, on average, younger stel- 
lar populations than pressure supported systems. This evidence 
suggests that rotationally supported dEs are low luminosity late- 
type galaxies which recently entered the cluster and lost their 
gas because of the interaction with the hostile environment, be- 
ing transformed, on short time scales, into dEs. It seems thus 
clear that not all dwarf early-type galaxies are the low luminos- 
ity extension of massive ellipticals. 

The new kinematic data in our hand support this scenario: 
rotationally supported systems have rotation curves similar to 
those of late-type galaxies of similar luminosity and follow the 
Tully-Fisher relation, the most representative scaling relation for 
late-type systems. Since the angular momentum of these objects 
is conserved, the most plausible scenario for gas stripping is the 
ram pressure exerted by the dense and hot IGM on the fragile 
ISM of the low luminosity star f orming galax i es freshly enter- 
ing the cluster environment (see Bos elli et al.l (|2008ab) for an 
extensive discussion). 

For these rotationally supported objects gravitational inter- 
action with the cluster potential or with other cluster members 
(galaxy harassment) can be excluded since they would, on rel- 
atively short time scales, reduce the angular momentum of the 
perturbed galaxies, leading to the formation of pressure sup- 
ported systems. This process, however, could still be invoked to 
explain the kinematic and structural properties of the remaining 
pressure supported dEs populating the core of the cluster (half of 
our sample), whose statistical analysis of their scaling relations 
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(Fundamental Plane) will be the subject of a future communica- 
tion. 

It is indeed possible that, as for the massive galaxies, gravita- 
tional interactions played a major role at early epochs, when the 
velocity dispersion of the cluster was lower since galaxies were 
accrete d through small groups (p reprocesing) than at the present 
epoch (iBoseUi & Gavazzill2006l) . 



9. Conclusions 

We present medium resolution (R~3800) spectroscopy for 21 
dwarf early-type galaxies, 18 located in the Virgo cluster and 3 
in the field. These spectra has been obtained with the IDS at the 
INT(2.5m) and with ISIS at the WHT(4.2m) at El Roque de los 
Muchachos Observatory (La Palma, Spain). We have used these 
data to measure the kinematic profiles of these systems and cal- 
culate their maximum rotation velocity as well as their central 
velocity dispersion. In order to guarantee the reliability of the 
data, we have compared our observations with the data of 972 
simulated galaxies in 5 different wavelength ranges correspond- 
ing to those used during the observations, and we have run, for 
each simulated galaxy, 100 Monte-Carlo simulations. The com- 
parison between observed data and simulations shows that the 
adopted data extraction technique is appropriate for measuring 
the kinematic parameters of the target galaxies. We have also 
shown that velocity dispersions can not be measured for S/N ra- 
tios below 15, while for radial velocities with S/N>10 accurate 
results are obtained. 

Our analysis has shown that dEs have on average dynami- 
cal M/L ratios within the effective radius smaller than those of 
massive ellipticals and dSphs (in average log(Ti) = 0.3 + 0.0 
log(Yi )). We thus confirm that, within the effective radius, dEs 
are not dark matter dominated objects. 

We have found that rotationally supported dEs have rotation 
curves similar to those of star forming systems of similar lu- 
minosity and follow the same Tully-Fisher relation. Combined 
with the evidence that these systems are young objects with 
disk- like structures gen erally located in the outskirts of the clus- 
ter ([Toloba et al. 2009), these observations are consistent with a 
picture where these rotationally supported dEs result from the 
transformation of star forming systems that recently entered the 
cluster and lost their gas through their interaction with the envi- 
ronment. The observed conservation of the angular momentum 
in the rotationally supported dEs suggests that a milder ram pres- 
sure stripping event as the responsible of the gas removal has 
to be preferred to more violent gravitational interactions (ha- 
rassment) which would rapidly heat up the perturbed systems. 
Therefore, all these evidences suggest that dEs are not the low 
luminosity end of massive early-types because if that was the 
case all dEs would be rotating with v max /cr higher than those of 
Es, but a population of non-rotators has also been found and, 
in addition, the evidences of being stripped late-type spirals are 
strong enough as to consider it as a possible origin of dEs in 
clusters. 
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Fig.A.l. Simulations performed to quantify the bias that can 
be introduced in the measurement of v max . Plotted are the 100 
folded rotation curves computed with random numbers dis- 
tributed as a Gaussian with width the typical errors of the target 
galaxies for that radius. Red squares and black dots are the left 
and right arms of the unfolded rotation curves respectively. The 
thick black line is the mean v max for all the simulations. The grey 
shaded area is the scatter for this mean. 
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Appendix A: Quantification of the bias introduced 

in v max 



Using the technique described in Section 15.21 to measure the 
maximal rotation speed, it seems that none of the galaxies have 
zero rotation (look at the v max presented in Table [4]). However, 
rotation curves like PGC 1154903 or VCC1087 (Figure |6]) ap- 
pear to be statistically consistent with no rotation. The fact that 
we always find positive maximum velocities is a consequence of 
the method used to measure it. To quantify the bias introduced 
by using this technique we have run some simulations. We have 
taken a zero-rotation object with errors typical of those galaxies 
that statistically are non-rotators. We have taken a typical rota- 
tion curve with 1 1 bins, at 0, 2", 4", 7", 12" and 16", with sym- 
metrical errors of 2, 5,7, 10, 12 and 15 km s _1 respectively. From 
these errors we have generated 100 simulated rotation curves as- 
signing to each radius a random number with a Gaussian distri- 
bution, of which the width is the error associated to that radius. 
After folding these simulated rotation curves, shown in Figure 
IA.1I we calculated v max following exactly the same technique as 
the one we used for the target galaxies. We then obtained a mean 
value for the 100 values of v max and its scatter, 9 ± 6 km s , 
shown as a thick black line and a grey shaded area in FigurefAT) 
As a result, we consider that those galaxies with v max < 9 km s _I 
are not rotating, based on our data. For the other galaxies the 
rotation is significant (three times above the standard deviation, 
except for VCC1122 and VCC1261, see Table H 3rd column), 
so that the systematic bias described here is much less relevant. 
For VCC1122 (17.3 ± 7.7 km s 1 ) and VCC1261 (13.9 + 5.2 km 
s _1 ) the rotation is marginal, as also shown by their low v max /cr. 



Appendix B: Absolute I band magnitudes, optical 
and half-light radii and ellipticity 

The /-band images for our sample of 21 d warf galaxies wer e 
drawn from Sloan D igital Sky Survey (SDSS.lYorket al.1 (12000)) 
data release 6 (DR6. lAdelman-McCarfhv et al.1 (120081) 1 and con- 
verted to Johnson-Cousins systems following Appendix ICl The 
photometric parameters were calculated using the IRAF task el- 
lipse. 

To remove the stars from the images we used the IRAF task 
fixpix. This task allows us to remove the stars interpolating the 
surrounding galaxy area. To improve the final outcome, we av- 
eraged the results of interpolating along the horizontal and ver- 
tical directions for each star. For those galaxies that were on 
the edge of the FITS images or had a very bright nearby star, a 
more careful procedure was implemented. Taking advantage of 
the fact that the galaxies are ellipticals, so that they have smooth 
and axisymmetric surface brightness profiles, the bmodel task of 
IRAF can be used to replace the affected areas of the galaxy by 
the azimuthal average of the unaffected ones. The output from 
bmodel was used only to replace a small fraction of pixels, so 
this procedure was not affected by the possible presence of more 
subtle features, such as bars or spiral arms, which could not be 
reproduced by the model provided by the bmodel task. Once ex- 
tremely bright nearby stars had been removed the same proce- 
dure as above was followed with ellipse and fixpix. 

The procedure followed to run ellipse is dependent on the 
parameters we want to measure. First of all we run ellipse fixing 
only the center of the galaxy assuming a step between isophotes 
of 1 pixel, the rest of the parameters were left free. We also made 
the masks for the stars to be removed as described above. With 
the aim of measuring the absolute magnitude, R opt and R e ff we 
run ellipse again fixing, the center of the galaxy, the ellipticity 
and the position angle (PA) to avoid overlap between consecu- 
tive isophotes. The adopted e and PA in this case are the typical 
values in the outer parts of the galaxy (beyond l.5-2R e ff, region 
where these two parameters stabilise). To measure e and C4 we 
run again ellipse after removing the stars leaving fixed only the 
center of the galaxy. 

Asymptotic magni tudes and the radii were derived as in 
iGil de Paz et al.1 (120071) . We first computed the accumulated flux 
and the gradient in the accumulated flux (i.e., the slope of the 
growth curve) at each radius, considering as radius the major- 
axis value provided by ellipse. After choosing an appropriate 
radial range, we performed a linear fit to the accumulated flux as 
a function of the slope of the growth curve. The asymptotic mag- 
nitude of the galaxy was the Y-intercept, or, equivalenly, the ex- 
trapolation of the growth curve to infinity. Once the asymptotic 
magnitude was known, the optical and effective radii of each 
galaxy were obtained as the major-axis of an elliptical isophote 
containing 83% and 50% of the total flux respectively. For the 
asymptotic magnitudes different sources of error have been con- 
sidered (see Appendix 0. The resulting uncertainty is ~0.02 
mag. 

The ellipticities (e) were measured as the mean value be- 
tween 3" and the R e ff, the galaxy region covered by our spectro- 
scopic observations. 



Appendix C: Errors in magnitudes 

The zero points (ZP) and the errors in the /-band magnitudes 
have been computed as described in SDSS documentation. The 
ZP have been obtained from F 0, the flux a source produces in 
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counts per second in the image, calculated as a function of three 
parameters (aa, kk and airmass) defined as: 



L exp 



j qO A(aa+kkxairmass) 



(Cl) 



where the exposure time (t exp ) is the same for all the SDSS im- 
ages (53.91 seconds). The uncertainties in the /-band magnitudes 
are affected by different sources of error: firstly, the errors in the 
flux, that can be calculated following the equation: 



AF 



F + sky 



gain 



+ Npix(dark variance + Asky) 



(C2) 



where F is the total flux in counts, the sky and Asky are the 
background sky and its error (in counts), the gain and the 
dark variance are given in the header and N p i x is the number 
of pixels in the largest aperture where the flux is measured. This 
error was typically 10~ 3 mag. Other error sources are the error 
introduced in the fit to the growth curve (between 10~ 3 mag and 
6 x 10~ 3 mag), and the error due to photometric zero point dif- 
ferences between the different scans of SDSS, which might lead 
to an error of 0.01 mag. SDSS /-band magnitudes are not ex- 
actly in the AB system, so an error of 0.01 mag might also 
be introduced (see SDSS documentation about the photomet- 
ric flux calibration). And finally, we have transformed our data 
from the SDSS /-band to t he Johnson-Cousins / -band assuming 
mi = mi - 0-52 + 0.01 mag Fukugi ta et alj (Q995) given that their 
r — i colour ranges from 0.23 mag to 0.57 mag. 

Adding quadratically all these sources of error, the final esti- 
mated error is 0.02 mag for the apparent /-band magnitudes. 



Appendix D: The C 4 Boxyness/Diskyness 
parameter 

The boxyness/diskyness parameter is defined as the fourth mo- 
ment in the Fourier series as follows 



/((D) = /„ + Y[S k sin(W) + C k cos(k<b)] 



(DA) 



/(<!>) is the intensity measured in each isophote. The first two 
moments in this series describe completely an ellipse. Higher 
order moments (k > 3) define deviations of the isophotes from el- 
lipses. The third order moments (S3 and C3) represent isophotes 
with three fold deviations from ellipses (e.g. egg-shaped or 
heart-shaped), while the fourth order moments (S 4 and C4) 
represent four fold deviations. Rhomboidal or diamond shaped 
isophotes have nonzero S4. For galaxies that are not distorted 
by interactions, C4 is the most meaningful moment indicat- 
ing the disky / boxy shapes of the isophotes (see Figure 1 from 
iPeletier et alj d 19901) for an example of these different shapes). 

C 4 is measured in /-band SDSS images using ellipse that 
performs equation lD. ll along the radius of the galaxy fixing only 
the center of the galaxy and leaving the rest of ellipse parameters 
free, as described in Section loTTI Figure lDTl shows C4 as a func- 
tion of radius for 3 dEs. Due to the large changes of C 4 with radii 
taking an averaged value is therefore not the best way to detect 
disks in these galaxies, especially if they cover only a limited 
range in radius. We have thus adopted the following procedure: 

If at least one prominent bump is detected, which has a width 
larger than ~6" inside three effective radii (above this radius the 
scatter of C4 and its error becomes too large as to be reliable), 
we consider the galaxy to be disky, and assign the maximum C4 




Fig. D.l. Examples of C 4 vs. radius for three galaxies: VCC397, 
VCC308 and VCC1695. The dashed purple lines indicate the 
region between 3" and the 3R e /f. The dotted purple line shows 
the R e ff- In the upper panel one can see that if an average value 
is used between the dashed lines, C4 will be compatible with 
zero and as a consequence, the prominent disky structure will be 
smeared out by the adjacent regions. In contrast, in the middle 
panel a galaxy with no clear disky structures is shown. In this 
latter case the errors are larger. The bump in C4 only covers ~ 
3", while the rest of the galaxy is boxy. In this situation a mean 
value of the C4 and its scatter is more representative. In the lower 
panel VCC1695 is an example of a typical boxy shaped galaxy. 



between 3" and 3R e /f to the global C 4 . The error in this mea- 
surement has been estimated by dividing the photometric error 
of the maximum of C\ by the square root of the number of points 
that describe the disky structure in order to quantify the reliabil- 
ity of the bump considered. If the bump is described by a large 
number of points it is highly likely that the bump is truly there 
and as a consequence the error will be small, but if the number 
of points is small but the photometry is of high quality then the 
error will be small again. In any other case the error will be large 
and the result must be used cautiously. 

Otherwise, if the values oscillate around C4 ~ (lower panel 
of Figure |D~TT ), are always negative or if there is a bump with a 
small radial coverage (see middle panel Figure |D~TT >, we assign 
a mean value and its scatter between 3" and 3R e ff to the global 
C4. In this case, the RMS quantifies simultaneously the quality 
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of the photometry and the possible presence of small bumps (as 
it is the case of VCC308, middle panel of Figure iDTTb . 

The results obtained for this parameter are listed in Table [6] 
and plotted vs. the anisotropic parameter, (v max /<r)*, in Figure 
[TOl The agreement between the C\ classification i n disky/boxy 
galaxie s and the morphological classification from lLisker et alj 
(2006fl) is evident but apart from three red dots. These filled 
circles correspond to VCC917, a rotationally supported galaxy 
(above the horizontal dash ed line) with strong disky isophotes 
but no structure found bv iLisker etafl (l2006ah ): VCC1122, a 
not rotationally supported dE but with an appreciable rotation in 
Figure[6]an d very important di sky structures in the inner R e jf not 
detected bv lLisker etaP (1200631) 1: and VCC1912, inside half the 
effective radius a moderate rotation is found in this system with 
a clear disky feature that peaks at \.5R e ff. For this galaxy, how - 
ever, no underlying structure was found bv lLisker et al.l (12006a). 
More importantly VCC308 and VCC856, two rotationally sup- 
ported galax i es with bo xy C4, present promine nt spiral arms in 
iLisker et alJ (l2006al) . In lFerrarese et al.l (120061) . based on ACS- 
HST images, VCC856 also shows spiral arms but their analysis 
of the isophotes' shapes shows that they are boxy too. Looking 
at Table [6] we see that both galaxies are nearly face-on, which 
means that the isophotes are boxy since face-on disks are round 
and not disky. As a consequence we emphasise the fact that boxy 
isophotes could miss disk features (mainly if the galaxies are 
face-on), but not the other way round (see VCC917). 
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Fig. 7. Comparison of our kinematic profiles with other works. In each diagram it is shown the rotation curve (upper panel) and 
the velocity dispersion profile (lower panel). The bottom X-axis is measured in arcseconds and the upper X-axis is measured as 
a fraction of the effective radius (R e ff) of each galaxy in i band (see Section 6). For Van Zee et al. (2004) we only present then- 
velocity profiles for Mgbl, more similar in wavelength to our data. 
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